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Abstract: We consider the adjacency matrix of the ensemble of Erdds-Renyi random graphs 
which consists of graphs on N vertices in which each edge occurs independently with probability 
p. We prove that in the regime pN ^ 1 these matrices exhibit bulk universality in the sense that 
both the averaged n-point correlation functions and distribution of a single eigenvalue gap coincide 
with those of the GOE. Our methods extend to a class of random matrices which includes sparse 
ensembles whose entries have different variances. 


1 Introduction 

The universality of the spectral statistics of random matrices has been a central subject since 
the pioneering works of Wigner [33], Gaudin [22], Mehta [26] and Dyson [7]. The hrst such 
universality result is the global semicircle law of Wigner which states that under some weak 
moment conditions, the empirical eigenvalue distribution of a matrix with i.i.d. entries converges 
weakly to the deterministic semicircle law 

1 ^ 1 ^_ 

^h\E\<2W^-E^ ( 1 . 1 ) 

in the appropriate scaling. 

The Wigner-Dyson-Gaudin-Mehta conjecture, or ‘bulk universality’ conjecture, states that 
the local statistics of the eigenvalues of random matrix ensembles should be universal in the sense 
that they depend only on the symmetry class of the random matrix ensemble but are otherwise 
independent of the law of the matrix entries. Here, local statistics refers to the behaviour of the 
eigenvalues in the scaling in which their typical distance is order 1. 

A prominent class of random matrices are Wigner matrices. These matrices have independent 
centered entries with a uniform subexponential decay condition and identical variances. The 
Wigner-Dyson-Gaudin-Mehta conjecture for Wigner matrices was recently established in a series 
of papers [9, 10, 12,14,16, 19] for all symmetry classes. Parallel results in various cases were 
obtained in [31,32], 

The conclusion of the papers [9,10,12,14,16,19] was that Wigner matrices, and even the wider 
class of generalized Wigner matrices in which the variances of the entries may differ, exhibit bulk 
universality in the following two forms. The first is that the n-point correlation functions are 
universal after averaging over a small energy window. The second is that the distribution of 
the eigenvalue gaps with a fixed label are universal. For Wigner matrices, the universality of 
the averaged n-point correlaton functions is equivalent to the universality of a local average of 
eigenvalue gaps. However, there is no rigorous mathematical relation between the universality 
of the eigenvalue gaps with a fixed label and the universality of the n-point correlation functions 
at a fixed energy. Universality at a fixed energy has recently been established for all symmetry 
classes in [5], but we will not be concerned with this type of convergence in this work. 

Wigner matrices were orginally introduced in [33] by Wigner to model the spectra of heavy 
atoms, and are widely used to model systems in which all elements strongly interact with one 
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another. However, for systems in which the links between different elements are broken, a better 
description is offered by the so-called sparse (or dilute) random matrices which have an average 
of pN nonzero elements per row, for p <C 1. 

Aside from theoretical physics models, sparse random matrices also arise in graph theory in 
the study of sparse random graphs. Perhaps the simplest example is the Erdos-Renyi ensemble 
which consists of a random graph on N vertices in which each edge is chosen independently 
with probability p. The adjacency matrix of this graph is called the Erdos-Renyi matrix. The 
Erdos-Renyi matrix has typically pN nonzero entries in each column and is sparse if p <C 1. As 
the matrix entries take values 0 or 1, the mean of the entries is not 0. Ignoring the nonzero 
mean for the moment, the Erdos-Renyi matrix can be viewed as a singular Wigner matrix, as the 
probability distribution of the matrix elements is highly concentrated around 0. The singular 
nature of this ensemble can be expressed by the fact that the k-th moment of a matrix entry is 
bounded by 

( 1 . 2 ) 

When p <C 1, this decay in k is much slower than in the case of Wigner matrices. 

It was conjectured in [20] that for sparse random matrices there exists a critical value pc > 1, 
such that for pN > pc, the bulk eigenvalues are strongly correlated and are characterized by 
GOE/GUE random matrix statistics; for pN < pc, the eigenvalues remain uncorrelated and 
follow Poisson statistics. This conjecture is supported by a wealth of numerical simulations 
[20,24] and a nonrigorous supersymmetric approach [21,27]. The best rigorous result in this 
direction was obtained in the works [9,11] and asserts that if 

pN > Ar2/3+'^, (1.3) 

then the averaged n-point correlation functions of the Erdos-Renyi ensemble coincide with the 
GOE. 

In the present work we prove that in the regime 

pN > N^, (1.4) 

the local statistics of the Erdos-Renyi ensemble exhibit bulk universality. In addition to proving 
that the averaged n-point correlation functions coincide with the GOE, we also prove the uni¬ 
versality of the eigenvalue gaps with a fixed label. To further place the present work in context 
we recall the three-step strategy developed in [9,10,12,14,16,19] for proving universality for 
Wigner matrices: 

(1) Establish a local semicircle law controlling the number of eigenvalues in windows of size 
log{Nf/N. 

(2) Analyze the local ergodicity of Dyson Brownian motion (DBM) to obtain universality for 
Wigner ensembles with a small Gaussian component. 

(3) A density argument comparing a general Wigner matrix to one with a small Gaussian 
component. 

For an overview of this three-step strategy and more details we refer the reader to [17]. The local 
semicircle law for sparse random matrices in the regime (1.4) was established in [11]. However, 
in [9] Steps (2) and (3) were only completed for sparse random matrices in the regime (1.3). 

The key input from Step (1) into Step (2) is a high-probability a-priori bound on the eigen¬ 
value locations which is a corollary of the strong semicircle law. In the case of Wigner matrices, 
this bound is optimal and it allows one to conclude that local equilibrium is reached by DBM 
in times t = N^/N. Sparse random matrices do not obey as strong a semicircle law and so the 
time to equilibrium found in the work [9] was much longer. Moreover, due to the slow decay of 
the third moment, the approximation in Step (3) is not as strong in the case of Wigner matrices 
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and so could not be used for the large times required by Step ( 2 ). These two factors led to the 
condition (1.3) of [9]. 

In the recent work [25], the optimal time of Dyson Brownian motion to local equilibrium 
was established for a wide class of initial data (see [ 8 ] for related results on DBM with general 
initial data). Using this as an input we will prove that DBM reaches local equilibrium in the 
optimal time t = N'"/N when the initial data is a sparse random matrix. For the comparison of 
correlation functions, Step (3) was obtained in [9] via a Green function comparison theorem. In 
this paper, we will use a lemma of [ 6 ] which asserts continuity of DBM when viewed as a matrix 
Ornstein-Uhlenbeck process. It is interesting to note that this continuity lemma provides a very 
convenient tool for Step (3) in the sparse setting whenever a “weak local semicircle law” is valid 
- and in the case of sparse random matrices this is provided by a result of [ 11 ]. 

The universality of a single gap was established in [16] for Wigner matrices, i.e., for p 0 ( 1 ). 
The work [25] also yields gap universality for DBM after the optimal time t = N’' ^ and so our 
task is similar to the proof of the correlation function universality in that we must establish 
Step (3) and compare the gap distributions. However, the completion of Step (3) presents a 
major difficulty. Previously, for gap universality this step was based on results of [31,32] or [23] 
which states that the gap distribution of two Wigner ensembles coincide provided that the first 
four moments of these two ensembles match. However, these results were based on two inputs: 
firstly, certain level repulsion estimates; secondly, an optimal eigenvalue rigidity estimate. 

Optimal eigenvalue rigidity estimates for Wigner ensembles were proven in [18,19]. This 
estimate states that for any eigenvalue Aj in the bulk we have that [A* — 7*1 < with 

overwhelming probability, where 7 * is the deterministic classical location of the i-th. eigenvalue. 
The best known rigidity result for sparse random matrices is from [11], where it was shown that 
the bulk eigenvalues satisfy [A* — 7 j| < with overwhelming probability. 

Moreover, we do not expect that optimal rigidity holds for sparse random matrices. In fact 
in [30], it was shown that for sparse random matrices the linear statistics 

N 

i=l 

converges to a normal random variable with variance 0{{N'^p)~^), for cj) satisfying some reg¬ 
ularity conditions. This implies that the fluctuations of the eigenvalues are at least of order 
on average, and so we do not expect optimal rigidity to hold if p <C 1 . 

As mentioned above, the lack of rigidity for sparse ensembles resulted in the longer time to 
equilibrium for DBM being found in [9], and it again causes difficulty in trying to compare gap 
statistics. Rigidity results are a crucial input in establishing level repulsion estimates for Wigner 
matrices which are needed in order to compare the gap statistics of two ensembles. It was proven 
in [25] that a level repulsion estimate will hold for DBM after a short time. We show that one 
can combine the delocalization of eigenvectors together with the Ornstein-Uhlenbeck continuity 
lemma of [ 6 ] to pass this level repulsion from DBM to the initial sparse random matrix. This 
level repulsion estimate then gives us a key input for Step (3) and we are able to conclude 
universality of the gap statistics. 

Previous level repulsion estimates were obtained in [5,13] for Wigner ensembles whose entries 
have a smooth distribution. Estimates without a smoothness condition were obtained in [31,32] 
and also in the very recent work [28]. A weak level repulsion estimate for Wigner matrices also 
follows from the results of [16]. 

In fact, our strategy outlined above applies to a wider class of random matrices than sparse 
or Wigner random matrices alone. We will prove that bulk universality holds for a class of 
random matrices obeying only a weak estimate on the distribution of its eigenvalues, a weak 
decay condition on the third moment of the entries and an eigenvector delocalization estimate. 

The remainder of the paper is outlined as follows. In Section 2 we introduce the random 
matrix models under consideration, which we will call ‘stable’ random matrices, and state our 
main results. In Section 3 we obtain bulk universality for Gaussian divisible ensembles. In 
Section 4 we state and prove our level repulsion results for stable random matrices. In Section 5 
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we complete Step (3) outlined above and compare the bulk statistics of a general stable random 
matrix and a Gaussian divisible ensemble. In Section 6 we prove that sparse random matrices 
are stable and conclude universality for sparse random matrices. 


2 Definition of model and main results 

In our paper we will only state and prove our results for real symmetric random matrix ensembles. 
All of our methods extend with only notational changes to complex Hermitian ensembles. 


2.1 Sparse random matrices 

In this section we introduce the class of sparse random matrices that we study. We follow the 
notations and definitions of [9, II]. The motivating example is the Erdds-R&yi matrix whose 
entries are independent up to the constraint that the matrix is symmetric, and equal to I with 
probability p and 0 with probability 1—p. It is notationally convenient to replace the parameter 
p with q dehned through 

q := ^/N^. (2.1) 


We allow q to depend on N. We also rescale the matrix so that the bulk of its spectrum lies 
in an interval of order 1. For the Erdos-Renyi matrix we dehne H to be the N x N symmetric 
matrix whose entries hij are independent up hij = hji and each element is distributed according 
to 


where we have defined 


7 f 1 with probability ^ 

9 1 0 with probability 1 — ^ ’ 



We further extract the mean of each entry and write 


( 2 . 2 ) 


(2.3) 


H = B + 'yq\e){e 


(2.4) 


where e is the unit vector 




(2.5) 


Note that the matrix elements of B are centered. It is easy to check that the matrix elements 
of B satisfy the moment bounds 

E[(.|] = i E||6yh<^^, k>2. (2.6) 

We are prompted to make the following definition. We introduce two parameters q and / which 
may be A^-dependent. 


Definition 2.1 (Sparse random matrices). H is a sparse random matrix with sparsity parameter 
q and mean / if it is of the form 


H = B + f\e){e\ (2.7) 

where / is a deterministic number satisfying 

0 < / < ( 2 . 8 ) 


4 



and B is a matrix with real and independent entries up to the symmetry constraint bij = bji 
which satisfy 


E[hj] = 0, E[|6,,f ] = 1 E[\bi,\^] < (2.9) 

for 1 < i < j < N and 2 < k < where C is a positive constant. We assume that 

q satisfies 

iV“ < g < ivV2 (2.10) 


for some a > 0. 


2.2 Universality of sparse random matrices 

Our main result is the bulk universality of sparse random matrices as defined above. 

Theorem 2.2. Let H he a sparse random matrix as defined in Definition 2.1, with sparsity 
parameter q satisfying 

<q< (2.11) 


for some number a > 0. Then H exhibits bulk universality in the following two forms. Firstly, 
H has the single gap universality in the bulk. For any k > 0 and index i € [[kN, (1 — k)A^]] 


lim E^^^[0{Npsci7i){Xi - Xi+i), ■ ■ ■ j dL Aj-|_n))] 

N^oo 

-E(^o^)[0(iVp,e(7*)(A*-Ai+i),--- ,iVp,,(7i)(Ai-A,+„))] =0. 


( 2 . 12 ) 


Secondly, the averaged n-point correlation functions of H are universal in the bulk. We denote 

(n) (n) 

the n-point correlation function functions of H and GOE by pfi and Pqq^ respectively, then 
for any (5 > 0 and E e (—2, 2), and b > 


f-E+h 
lim / 
N^ooJE-b 


I O(oi, • • • , On 

JR^ 

^ f pi I 


Psc{E) 


Psc{Ey 

CXl 

Npsc{E) 


-/h if + 


Ol 


r--E' + 


Npsc{E) 




Otr. 


Npsc{E) 


Npsc{E) 


dai...dan—= 0. 
2b 


(2.13) 


where the test function O G 


oo(]^n) 


2.3 Stable random matrices 

While our main goal is to study the sparse random matrices defined above, we note that our 
analysis applies to a somewhat more general class of random matrices. We consider an x 
real symmetric random matrix H = (/iij)i<i jxat, which satisfy E[/ijj] = / and E[{hij — f)'^] = Sij. 
We assume that there are constants ci and C 2 such that 

ciA^“^ < Sij < C2N~^, (2.14) 

and the mean 0 < / < may depend on N. 

We define the following matrix stochastic differential equation [6] which is an Ornstein- 
Uhlenbeck version of the Dyson Brownian motion. The dynamics of the matrix entries are given 
by the stochastic differential equations 


- f) 


d.Bij{t) 

Vn 


1 

2N Sij 


{hij{t) - f)dt. 


(2.15) 
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where B is symmetric with {Bij(t))i<i<j<]\f a family of independent Brownian motions. We 
denote Ht = {hij(t))i<ij<N, and so Hq = H is our original matrix. More explicitly, for the 
entries of Ht, we have 

hij{t) = f + e {hij{0) - f) + J dBij{s). (2.16) 

Clearly, for any t > 0, we have K[hij(t)] = /, and E[(hij(t) —/)^] = stj. More importantly, 
the law of hij{t) is Gaussian divisible, i.e. it contains a copy of Gaussian random variable with 
variance 0{tN~^). Therefore Ht can be written as 


where r = minj<j{A^Sjj}, G denotes a standard gaussian orthogonal ensemble, which is inde¬ 
pendent of h'^^ . The entries of the matrix h'^^ is given by 


(h[^A =/ + e (h,,(0)-/) + 

\ / ij 



1 + ^ij 


r(l - e r) 


t Bij{t) 


y/N 


where B is symmetric with a family of independent Brownian motions. 

We define the deformed matrix 6°‘^Ht by 




(2.18) 


where 9'^^ 
0 < 9f, = 


= 1 unless {k,l} 

< < 1 - 


{a, 6} in which case 0“^ = will be a number satisfying 


Definition 2.3. Let ^ be an x deterministic real symmetric matrix. We denote the 
eigenvalues of A as {Ai, A 2 , • • • , Aat} and corresponding eigenvectors {ui,U 2 , ■ ■ ■ , un}. For any 
(small) number <5 > 0, we call the matrix M (i-general if: 

(1) The eigenvectors of A are completely delocalized: supjj \ui{j)\^ < . 

(2) The eigenvalues of A do not accumulate: there is an universal constant C, such that for 

any interval I with length |/| > we have : Aj € /} < C|/|A^. 


Definition 2.4. We call the random matrix H stable if: 


(1) The entries of H are independent up to symmetry. 

(2) For any time t = where e > 0 can be arbitrarily small, the random matrix H^'^ 

defined in (2.17) satisfies the weak local semicircle law, i.e. for any (large) number D > 0, 
and (small) number <5 > 0, the following holds with probability larger than 1 — N ^, 


— TT{Hj:^^ — E — irj) ^ — msc{E -|- irj) 




(2.19) 


uniformly for —5 < E < 5, N < ?? < 10, where c{S) > 0 is some constant depending 
on 6, and spec(i7^^^^) C [—3,N^) for some fixed C. 

(3) There exists some universal constant a, such that — f\^] < 

(4) For any time 0 < s < t, any (large) number D > 0, and (small) number <5 > 0, 9°‘^Hs 
is (i-general with probability larger than 1 — N~^, with the constants in Definition 2.3 
uniformly in s. 

Remark 2.5. The second condition (2) implies that for any k > 0, the eigenvalues Xi{Hj:^'^) £ 
(—2 -|- 2 — K^/^) for i £ [[2kA^, (1 — 2«:)A^]] with overwhelming probability. 
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Remark 2.6. In order to simplify our proof we have assumed that the matrix elements are 
independent. Independence is mainly used in the comparison Lemma 4.3, which will still hold 
if the matrix entries of H are weakly correlated. 

Remark 2.7. The motivating example of our paper is the sparse random matrix, and we have 
therefore assumed that the Stieltjes transform of the empirical eigenvalue distribution of 
is close to nisc, i-e., the semicircle law. The semicircle law, however, does not play an active 
role and our methods can be applied to the case in which the semicircle law is replaced by other 
densities. We will not pursue this direction and refer the interested reader to [2, 3] and [1] for 
examples in which the limiting eigenvalue density differs from the semicircle law. 

In this paper we will prove that the local statistics of stable random matrices are universal. 


Theorem 2.8. Let H be a stable random matrix as defined in Definition 2.4- The local statistics 
of H in the bulk are universal. Firstly, H has gap universality with a fixed label in the bulk. For 
any k > 0 and index i G [[kN, (1 — «:)A^]], we have 


lim E^^'>[0{Npsci'yi){Xi - Ai+i ),••• , N Aj-i-n,))] 

vf —>-00 


The averaged n-point correlation functions of FI are universal in the bulk. We denote the n-point 

fn) (n) 

correlation function functions of H and GOE by pfi and Pqq^ respectively, then for any d > 0 
and E G (—2, 2), and b > 




Qfl 


[EYPgoe + 


Otl 


PsciEfi 

Above, the observable O G 


Npsc{E) 


,---E' + 


NpsciE) 


Oir, 


,---E' + 


Oir, 


NpsciE) 


Npsc{E) 


dai...da„^ = 0. 
2b 


( 2 . 21 ) 


3 Bulk universality of Ht 

The goal of this section is to establish bulk universality for the matrix valued stochastic process 
Ht defined as in (2.17) after a short time t = 

Theorem 3.1. Let H be a stable random matrix, and let Ht be defined as in (2.17). For any 
small e, K > 0 there is a constant c > 0, which depends on e, such that the following holds for 
t = and any index i G [[kN, (1 — k)A^]] 


E(^‘)[0(iV(Ai - Xi+i), • • • , N{Xi - Ai+„))] 


< 


Moreover, for any 5, k > t), E G (—2 + k,2 — k) and b > N we have 


rE-\-b 

lE-b 

1 


PsciEf 


/ 0(ai,...,c 

JR^ 

:,Pgoe (^E' + 


P.c{Ey 


Pm ( + 


ai 


NpsciE) 


,...,E' + 




NpsciE) 


Ol 


NpsciE) 


,...,E' + 


Olr. 


Np^E) 


dQ;i...da, 


dE', 


< N- 


where the test function O G C'“(M”). 


(3.1) 


(3.2) 
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For the proof we shall first restate the main result of [25] in a form convenient for our proof. 
For this we will introduce some notation. For a deterministic matrix A we define 

I r(\ — 1 

At:=A + ^tG, ^^ (3.3) 

where G is a GOE matrix and r is the constant from (2.17). We denote by mt the Stieltjes 
transform of the free convolution of the empirical eigenvalue distribution of A and the semicircle 
law of and so mo is the Stieltjes transform of empirical eigenvalue distribution of A. More 
explicitly, mt is defined as the unique solution to the functional equation 

1 ^ 1 

mt{z) := mo{z + i^jmtiz)) = ^ >0,lm[z] >0. (3.4) 

The free convolution is well-studied, see, e.g, [4]. It is known that mt is the Stieltjes transform of 
a measure with a density which we denote by pt which is analytic on the interior of its support, 
for any t > 0. Denote the classical eigenvalue locations of the density psc and pt by 7 * and 
respectively, 

/ 7i ^ ni,t j 

Psc{x)dx = —, J Pt{x)dx = —. 

The following follows from Theorem 2.5 of [25]. 

Theorem 3.2. Fix e > 0 and suppose that there are constants ci > 0 and Ci > 0 such that 

Cl < Im[mo(E + iy)] < Ci (3.5) 

for all E G (—2+k, 2 —k) and < 7 < 10. Suppose furthermore thatspec{A) C [—C,N^) 

for some fixed C. Let i be such that A* (A) G (— 2 + 2«;, 2 — 2k). Then for t = there exists 

a small constant c> 0, which depends on e, such that for indices ii, ...,in £ o,nd i^ < N'^, 

|E(^‘)[ 0 (pi( 7 ,,i)iV(A, - A,+,J, • • • ,ptht,t)NiX, - Ai+,J)] 

,^^^( 7 .)JV(A, _ A,+,J)]| < (3.6) 

for all large enough N > N{e). Above j is any index satisfying j G [[kiA^, (1 — ki)A^]] where 

Kl > 0 . 

We will apply the above theorem by conditioning on and taking A = h[^\ However, 
we must take care of the scaling by pti'yi,t) in the above theorem statement. This becomes a 
random variable depending on We want to replace it by the deterministic quantity Psci'Ji)', 
we will see that our assumption on the weak local law of Flj:^'^ combined with continuity of the 
free convolution allows us to do this. 

Define the spectral domain 

F = {z = E + ip:Ee (-5, 5 ), < ^ < 10 }, (3.7) 

and the set of real symmetric N x N matrices 

V := {A : lmo(2:) — mgc(z)| < N~‘^, z G F} Pi {spec(A) C (—3, A^^)}. (3.8) 

where oj satisfies 

0 < a; < min{e/3, c(e/3), (1 — e)/4}, (3.9) 

and c(e/3) is from the constant in (2.19) of Definition 2.4. 

Lemma 3.3. Let k,u] > 0 and A^V as above. Then, for any e > 0 and t = 

\mt{z) - msc{z)\ < CN~‘^, (3.10) 

uniformly for any z E {z = E + ip : E ^ (~4,4), A^“^+^/3 < 77 < 9}. 
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Proof. From Jensen’s inequality we have: 


\mt\ S 


< — 






N ^ |Aj(y4) — z — lm[z + '&‘fmt{z)] 


<^t 


-2 


Therefore we have < 'dt = Since A G P, for any z = E+irj, such that E G (—4,4) 

and < rj < 9, we have that 2 + 'd^mt G E. The defining relation (3.8) leads to 


mt{z) =mo{z + 'd^mt{z)) =msc{z + '&t'^tiz))+ 0{N ‘^) (3.11) 


To control msc{z + we have the following stability estimate of rUsc'- for any z, lS.z with 

non-negative imaginary part, and |Azl < 1, then 

\'<^sc{z + Az) — msc{z)\ < 2\/S.z\^^’^. (3.12) 


Therefore (3.11) gives us 

\mt{z) - msc{z)\ < \mt{z) - msc{z -h 'djmt)\ + \msc{z + rl^mt) - msc{z)\ < CN~‘^, 
given that w < (1 — e)/4. □ 

From this lemma we conclude the following. 

Lemma 3.4. For any e,K > 0, time t = any real symmetric matrix A £ V, and 

E G (—2 + k,2 — k) we have 


\pt{E)-psc{E)\<CN-‘^. (3.13) 

Moreover for any index i such that Aj(A) G (—2 + k,2 — k), we have 

\pt{Et)-pM\<CN-^/\ (3.14) 

where the constant oj is from the definition (3.8) of the set T>. 

Proof. From [25, Lemma 7.1], the derivative of pt satishes \p'^{E)\ < C/t for E G {—2 + k,2 — k), 
where we use the fact = Oft). Moreover the same lemma also shows that pt{E) < C on 
(—2 -|- k/2, 2 — k/2). Given E G (—2 -|- k, 2 — n), we denote the intervals Ii = [E — E + 

A-1+2^/3], l 2 = {x-. iV-i+2^/3 <\x-E\ < k/2} and Jg = M\(/i U h). We take p = 
then 


\Pt{E) 



-Im [mt{E + ip)]\ 

TT 


{pt{E) - pt{x))pdx 
{x - E)^ + 7/2 


+ 


]_ f pt{E)pdx 

^ Jhuh (a: - P)2 + 7/2 


pt{x)pdx 
{x — P)2 -I- 7/2 


1 r pt{x)pdx 

TT ijg {X - F;)2 7/2 


< sup \pt{x) - ptiE)\ + ClV-^/3 ^ + c^ [ Pt{x)dx < CN-^/^, 

X£li ^ J/3 


(3.15) 


where we have used \pt{x)\ < C on Ii U I 2 . Moreover, we have 

-Im [mt{E + 77/)] - Psc{E) 
vr 

<— jim [mt{E + 77/)] — Im [msc{E + 77/)]] H—jim [msc{E + ip)] -Im [777sc(F')]l 

TT TT TT 

<N-‘^ + -yEj < CN-‘^, (3.16) 

TT 

given cj < (1 — e)/2. (3.15) and (3.16) together lead to (3.13). 
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Moreover, by our hypothesis that A £ D, spec( 74 ) is bounded below by 3. Therefore the 
density pt is also bounded below: supp pt C [—7/2,oo). Therefore if i is such that Ai(^) G 
(—2 + K, 2 — k) then 

\li,t-li\<N-^/\ (3.17) 

This follows from [25, Lemma 7.17], using Lemma 3.3 as input. Therefore, 

- Psc(7i)l < \pt{li,t) - pschi,t)\ + Ipscilpt) - Psc(7i)l < (3.18) 

□ 


3.1 Proof of Theorem 3.1 

By the hypotheses of stability of H we have that G V with probability greater than 1 — N~^ 
for any large D. We denote the density of free convolution of and 'PtG as p^, and its f-th 
classical eigenvalue location as 7 /^. From Theorem 3.2 we have 

(7M)iV(A* - A,+,,), • • • ,pf (7,^)iV(A, - 

-E(^o^)[C)(psc(7i)A^(A* - Ai+iJ,--- ,psc(70^(Ai - Ai+iJ)]| < iVT (3.19) 

But then by Lemma 3.4 we have that {p^i'Jijt) ~ Psci'yi)\ < therefore 

(7M)iV(A. - A,+,,), • • • ^pHl^,t)N{h - 

-E(^‘)[0(psc(7i)A'(Ai - A,+iJ,--- ,psc(7i)A'(Ai - Ai+iJ)| 77 f)]| < CN-'^G (3.20) 


for some C > 0 depending on first derivative and the support of the test function O. We 
therefore obtain that for G P that 

|e(^‘)[O(psc( 70A'(A* - Ai+,J,--- ,Psc( 70^(A* - Ai+,J)|i7f)] 
_E(GOi?)[o(p^^(^.)iV(A, _ ,p,,(7i)iV(Ai - Ai+iJ)]| < (3.21) 


if we choose c small enough, such that c < ui/2. And if we take expectation over h[^\ (3.1) 
follows. 

Now let be the re-point correlation functions of Ht conditioned on It is well known 

that the estimate (3.19) together with the optimal rigidity of the eigenvalues of Ht, [25, Theorem 
3.3] implies that for gV we have 


rE-\-b 
JE-b . 


0 ( 0 ^ 1 , Q^n) 


(n) 


Psc{Ey 


Psc{E) 


^ „(”) (r' 4- 

'^yPcOE + 


Oil 


iVpsc(^)’' 


:Pi 


.,S' + 


E' + 


ai 


NpsciE) 


,...,E' + 




Ckr. 


Np^E) 


dQ;i...dar 


NpsciE) 

dE' 


< N- 


(3.22) 


for b = N for any 5 > 0. For this argument, we refer the reader to, e.g., [14, Theorem 2.1]. 
We conclude (3.2) by integrating over Hj:^K □ 


4 Level repulsion for stable random matrices 

In this section we prove the following level repulsion estimate for stable random matrices. It 
will be used for the comparison of the single gap statistics between H and Ht in Section 5. 

Theorem 4.1. Let H he a stable random matrix as defined in Seetion 2.3. Given any 0 < r < 
a/ 8 , any (small) number re > 0, and any index i G [[kIV, (1 — re)A^]], we have 

¥i\XiiH) - Xi+iiH)\ < V-1-) < N-G\ (4.1) 
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Remark 4.2. The above estimate suffices for the comparison of the single gap statistics of H 
and Ht- We have not tried to optimize the exponent —t/ 2, which is far from optimal. The proof 
below is easily modified to give — r + 5 for any <5 > 0. 

It was proven in [25] that a level repulsion estimate holds for the matrix Ht for t = 

To obtain the level repulsion estimate for H, we need to prove that the change of eigenvalues 
up to time t = is negligible. For this we will repeatedly use the following lemma which 

asserts continuity of DBM when viewed as a matrix Ornstein-Uhlenbeck process. It is a minor 
modification of [ 6 , Lemma A.2]. 

Lemma 4.3. Let H be an N x N real symmetric random matrix H = {hij)i<ij<N, where the 
hij’s are independent up to symmetry constraint hij = hji. Suppose that for constants ci,C 2 ,C 3 
its entries satisfy ^[hij] = f and IE[(hjj — /)^] = Sij, where ciN~^ < stj < C 2 N~^ and / > 0 
may depend on N. 

Define Ht as in (2.15). Let F be a smooth function on the space of real symmetric matrices 
satisfying 


sup sup E 

0<s<t,a<b 


{N^\hM - ff + N\hab{s) - 


< B 


Above, the deformed matrix 6°-^Hs is defined by {9°'’^H)ki = / + {hj-i — f), where = 1 

unless {k, 1} = {a, b} and 9^ = 9^j^ is a number satisfying 0 < 9f^ = 9^^ < 1. Then 


|E[F(Rt)] - E[F(Ro)]| <CtB. 


(4.2) 


Given a real symmetric matrix A, we denote its eigenvalues by {Ai, A2, • • • , Aat}, and corre¬ 
sponding eigenvectors {ui,U 2 , ■ ■ ■ ,U]\f}- If A* is a simple eigenvalue of A, we define Pi to be the 
orthogonal projection to the one-dimensional eigenspace corresponding to Aj, and the resolvent 
Ri{A) the unique real symmetric matrix inverting Xt — A on the range of I — PfiA), and vanishing 
on the range of PfiA). RfiA) can be written explicitly as 


RfiA) = Y, 






Moreover Ri{A) can be written as the following contour integral, 


R^{A) 


1 

27ri 





dz, 


(4.3) 


where we pick the contour tv to enclose only Aj. From the above formula it is clear that RfiA) 
is a smooth function on a neighbourhood of A if A, is a single eigenvalue. We refer the reader 
to the book [29, Chapter XII] for related properties. 

If Xi is a single eigenvalue of A, we define the quantity 

= = (4.4) 


This quantity plays an important role in [31], where it was observed that it captures quantita¬ 
tively the derivatives of A*(A). Here we write it in terms of the Green function, and prove it 
is stable under the DBM (2.15). As a result, based on the idea of Green function comparison, 
it can be used to derive the weak level repulsion estimate Theorem 4.1, once we know such an 
estimate for larger times, see Theorem 4.4. 

Since QfiA) is not well-defined on the space of real symmetric matrices (it will blow up when 
Xi{A) is not a single eigenvalue), we have to compose it with a cutoff function xm^ where M := 
N‘^'^ and the (small) constant r > 0 will be chosen later. We choose the cutoff function xm{x) 
which satisfies the following two properties: ( 1 ) xm is smooth, and the first three derivatives are 
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bounded by some constant C, i.e. |Xm(®)I) IXm '(®)|5 (^) interval [0,M], 

Ixm^x) — xj < 1, and for x > M, xm{x) = M. 

If Xi is a single eigenvalue of A, then in a neighborhood of A, XM{Qi{A)) is smooth; if A* is 
not a single eigenvalue of A, then in a neighborhood of A, XM{Qi{A)) is constant, which is also 
smooth. Therefore XM{Qi{A)) is a well defined smooth function on the space of real symmetric 
matrices. 

Our proof of Theorem 4.1 consists of three steps: 

1. We use the level repulsion estimate for Ht to conclude the estimate K[xM{Qi{Ht))] < 
CN^r/2 ^ ^ Ar-i+^ 

2. We compare IKlxMiQiiHt))] for t = and KlxMiQiiHo))] using the continuity 

Lemma 4.3. Since e > 0 can be chosen arbitrarily small we will obtain 

\E[xM{Qi{Ht))] - E[xM{Qi{Hom < CN\ (4.5) 


3. Theorem 4.1 then immediately follows from E[xM{Qi{Ho))] < and the Markov 

inequality. 

In the remainder of this section we denote the eigenvalues of Ht by Xi{t) < A 2 (t) < • • • < Aw(t), 
and so the eigenvalues of H are Ai(0) < A2(0) < • • • < Aw(0). 

The following level repulsion estimate is an immediate consequence of [25, Theorem 3.6]. 

Theorem 4.4. Let H be a stable random matrix as defined in Section 2.3. Given any (small) 
number w > 0, and (large) number D > 0, we have that 

P {\Xiit) - Ai+i(t)| < ON-^) < + N-^ (4.6) 

for any i G [[kN, (1 — «:)A^]] and t > 

From this we derive the following estimate. 

Lemma 4.5. Let H be a stable random matrix. Then, 

E[xM{Qi{Ht))] < (4.7) 

for t = for any e > 0 and all small r > 0. 

Proof. By our assumption with probability larger than 1 — N~^, the matrix Ht is d-general in 
the sense of Definition 2.3. Combining this with Theorem 4.4, we have 

P [Ht is A-general and \Xi[t) — Ai±i(t)| > 6N~^) > 1 — — N~^. 


Define a dyadic decomposition: 

Utt = {j^i-. |A,(t) - A,(t)| < Uoo = {j : < |A,(t) - Ai(t)|} 

Un = {j : 2-iiV-i+^ < |A,(t) - Xi{t)\ < 2-Ar-i+'5}, l<n< [logs ■ 

For those Ht which are A-general in the sense of Dehnition 2.3, |tAi| < C2^N^, for 0 < n < 
[logs -^1 • the event that Ht is d-general and |Ai(t) — Ai±i(t)| > 9N~^, we derive the estimate. 


riog2 


Qi{Ht) - ^s X] 


1 




+ < 3CW6I 


5/1—2 


(4.8) 


Therefore we have 


W E 1 < + v- 
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where oj > 0 can be any small number. Therefore 


E[xM{Q^{Ht))] < 1 + 3CN^e-^ + + N-^). 


If we take 6 = N then we get E[xM{Qi{Ht))] < if we take <5 < r/2, a; < r/4 and 

D large. □ 


This finishes the proof of the hrst step. In order to apply Lemma 4.3 for the second step, 
we need to control the third derivative of XM{Qi{(^°‘^Hs)). 

Proposition 4.6. Let A he an N x N deterministic real symmetric matrix. If A is 5-general 
in the sense of Definition 2.3 and 


Q^{A) 


1 


E 


1 

\\fiA)-\fiA)\^ 


<M = N'^^, 


(4.9) 


then 


dabQM) < A; = 1, 2,3. (4.10) 

for some constant C. 

Proof. We denote G = {A — z)~^ the resolvent of A and by Xj and Uj the eigenvalues and 
eigenvectors of A. Notice that (4.9) implies |Aj — Aj±i| > . The same dyadic argument 

leading to (4.8) yields 


E 


1 

I 


EE 

n>0jGUn 


1 

I 


< CiV^+^+^. 


(4.11) 


Also we have the trivial bound for higher moments 

\ k/2 

,, \ „ I k>2. (4.12) 

jAi-A^l^y 

We denote by V the matrix whose matrix elements are zero everywhere except at the (a, b) and 
( 6 , a) position, where it equals one. From the formula (4.3), 

= (4.13) 

Notice that = {—lYk\{GV)^G. By the Leibniz rule, (4.13) can be written as a sum of 

terms in the following form 

E Tr ^ (GV)>^^G8lXj^^dz f (4.14) 

where /ci + A ;2 + ^3 + A :4 = k. We will only prove (4.10) for k = 3] that is, we will prove 

\d^^Qi{A)\ < CiV5r+85 (4^15) 

The computations for A: = 1,2 are much easier. To evaluate (4.14), we need to compute the hrst 
three derivatives of A,(A) with respect to the (a, b)-th. entry of A. We use the following formula 
to compute the derivatives of Aj, 

Aj =-Tr (f — - dz 

2-ki J A-z 


E 


|Ai - A, 


ifc — 
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where the contour encloses only Aj. The k-th derivative with respect to (a, 6 )-th entry is 

za'^^Gdz = tllplTr f z{aV)'-adz. 

For k = 1, 2, 3 respectively we have 

N 


I ^ f zuWuj 

dabXi = (X 


o( 2 ) \ . _ 


27ri ^ J iXd — z) 

j=i d 

-2 ^ f z{uj^Vu*^f 




2711 


6 


} 


2E 


{v^Vujf 


J 1 J 2 - r.j^i 


27ri 




Xi-Xd ’ 

z (rr ji F ) {uj^ V u *^) V u*^) 

(Aji — z)^(Aj 2 ~ ^)(^i3 “ 


(4.16) 

(4.17) 


^ (^i rt *T^ \ {u^VUj){u^VUi 

(a,-a,,)(a,-a, 2 ) 






(A. - A,)2 


(4.18) 


Some straightforward but tedious integration similar to (4.16), (4.17),(4.18) reveals that 
f 3l 

Qi{A) is a sum of the following terms (for simplicity of notation we write Vjk ■= u*Vuk)'- 


E 

E 

31^32- 

31^2^'^ 


n 


o — V- v3v- ■ F F A V- 

— rjj, Vjj Vjj, V^ 


2 t/3 


(Ai-Aj)5^ 




(A,-A,J4(Y-A,2) 


JJ’ JJ' 


Q — T/- T/.2 . T/. . T/.2 W. . W.2 , T/..w2 


(A- -AjJ3(A._Aj 2)2’ ““ ’J2J2-'iii2, iji, -'n''iii2’ ''JUi-'ilia- 


E 

jl,J 2 ' 

E 


0 


^'li 2 ^'2i3 ^'3il 


n — V ■ F? ■ F F^ F F-^ ■ F ■ F^ 

. E 


^'li2 ^213 ^jzh 


Ji,d2G3- Aj^)^(Ai Xj^){Xi Ajg) (Ai Aj22(Ai Aj2)^(Aj Ajg) 

h <32 •337^'^ h <32 •337^'^ 

Indeed, if one takes a close look at the expression (4.14), the only singularity enclosed by our 
contour is Aj. Therefore by Cauchy’s formula, the integral is sum of terms with denominators: 
02 A. — Xj), as appearing in the above expressions. 

Since the eigenvectors of A are completely delocalized we see that 


\Vjk\ < 

This, together with the bounds (4.11) and (4.12) yields (4.15) 


(4.19) 

□ 


4.1 Proof of Theorem 4.1 

By explicit computation, we have 

[xM{Q^{0'^’’H,))) = x'M{dabQ3f + ^XM^bQAbQi + XM€bQ^■ (4-20) 

Notice that the right hand side of (4.20) vanishes unless Qi{6°‘^Hs) < M. Since H is stable, 
9°‘^Hs satishes all the assumptions in Lemma 4.10, with probability larger than 1 — N~^ for any 
large number D. Therefore 


(xMiQiiO'^^Hs))) I < CiVl25+9r^ > ^ _ jy-D_ 
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On the complement of the above event we will use the following deterministic bound 






We want to apply Lemma 4.3 with F = xm ° Qi- For this choice of F we can take B to satisfy 


B < sup E 

0<s<£ 

l<i,j<N 


{N^\hij{s) - ff + N\hij{s) - f\)CN 


12(5+9t 




where the last factor “ is from the third moment of hij{s). We conclude that 

\nXM{Qi{Ht))] - E[xM{Qi{Hom < CtiVl+125+9r-« ^ ^ ^e+12S+9r-a _ ( 4 ^ 21 ) 

Since the two numbers e,6 > 0 can be arbitrarily small, therefore we can choose r = {a — e — 
125)/S and (4.21) simplihes to 

\E[xM{Q^{Ht))] - E[xM{Qi{Ho))]\ < CNF (4.22) 

Hence, 

E[xM{Q^{Ho))] < (4.23) 

and the proof is easily concluded by the Markov inequality. □ 


5 Bulk universality of stable random matrices 

In this section we prove bulk universality for stable random matrices H, i.e.. Theorem 2.8 by 
comparing the local statistics between H and Ht- This will yield the theorem as Theorem 3.1 
shows that the latter ensemble exhibits bulk universality. In the following of this section we 
denote the eigenvalues of Ht by Xi{t) < X 2 {t) < ■ ■ ■ < AAr(t). And so the eigenvalues of H are 
Ai(0)<A2(0)<---<Ajv(0). 

We will obtain the gap universality of H from the more general comparison result below. 

Lemma 5.1. Let H he a stable random matrix and let Ht he defined as in (2.15). Take t = 
Ar-i+e_ Then for e > 0 small enough we have 

hm E[O(iVA,,(0),--- ,iVAi„(0))]-E[O(iVAq(t),--- ,AA,Jt))]=0, (5.1) 

N^oo 

for any ii, Z 2 , • • • fin ^ (1 ~ ^)A^]] bounded test function O G C'°°(M"') whose first three 

derivatives are bounded. 

Proof. For simplicity of notation, we only state the proof for n = 1 case, i.e. for any i G 
[[kV, (1 - k)N]] 

lim E[O(VA*(0))] - E[0(iVAi(t))] = 0. (5.2) 

N—too 

Take a cutoff function pM such that pm{x) = 1 for x < M and pm{x) = 0 for x > 2M, 
where M = and r > 0 is a small constant. By the level repulsion of H and Ht from the 
previous section, we know that 

EiQfiHs) > N‘^^) < s = 0fi. 

Since O G C'°°(M) is bounded, we have that 

\E[O{NXfi0))]-E[OiNXi{t))]\ 

< \E[0{NXfiO))pM{Qi{Ho))] - E[OiNXiit))pM{Q^iHt))]\ 

+ ||0||oo {E{Q^{Ho) > N^fi+E{Qi{Ht) > 

< \E[0{NXfiO))pMmHo))]-nO{NXfit))pM{Q^{Ht))]\ + 
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Notice that 0{NXi{A))pM{QiiA)) is a well defined smooth function on the space of symmetric 
functions. Moreover, if the matrix A is h-general in the sense of Definition 2.3 the same argument 
as in Proposition 4.6 implies 


d^J}0{NXi{A))pM{Qi{A)) < 


where c and C are constants. Therefore by Lemma 4.3, we have 

\E[0{NXi{0))pMiQi{Ho))] - E[0{NX,it))pM{Q^{Htm < CiV'=('5+D+^-« ^ q, 


if we take c{S + r) + e < a/2. □ 

For the universality of the correlation functions, due to the lack of optimal rigidity, one 
can not directly deduce universality from (5.1). We need to prove the following Green function 
comparison lemma: 

Lemma 5.2. Let H be a stable random matrix as defined in Section 2.3, and let Ht be defined 
as in (2.15). Let 6 > 0 be arbitrary and choose an p with N~^~^ < p < N~^. For any 
sequence of positive integers ki, k 2 ,. ■., kn, any set of complex parameters zip = Ep ± p, where 
1 < J < km, 1 < m < n, \Ep\ < 2 — k, and the ± signs are arbitrary, we have the following. 
Let Gt{z) = {Ht — z)~^ be the resolvent and let F{xi,X 2 , ■ ■ ■ ,Xn) be a test function such that 
for any multi-index a = (ai, • • • ,an) with 1 < [a] < 4 and for any w > 0 sufficiently small, we 
have 


max 


|9“T(xi,X 2, • • • ,x„)| 


max \xj\ < N‘^ 




(5.3) 


and 


max 


\d°‘F{xi,X 2 , • • • , Xn)\ : max \xj\ < > < N^° 


(5.4) 


for some constant Cq. Then the following holds: 


ki 

E[F ( iV-'^iTr 




—kt 


Tr l[Gt{zl/) 1 -E[F{Gt^Go)]] 
i=i 


< (5.5) 


where c and C are constants depending on k, n, ki,k 2 , - ■ ■ ,kn and Cq. The second term above 
is the same as the first term, but with Gq replacing Gt everywhere. 

Proof. For simplicity of notation, we state the proof only for n = 1 and fci = 1 case, i.e. 

\E[F (iV-iT\'Gt(z))] -E[F(iV-iTrGo(z))]| < 


We will prove this lemma using Lemma 4.3. We must compute derivatives of the trace 
of the Green’s function of the deformed matrix We denote the resolvent of O^'^Hg by 

G{z) = {9°‘^Hg — z)~^. For the derivatives of G{z), we have | = |(—l)^/c! Tr(GP)^G|, 

where V is the matrix whose matrix elements are zero everywhere except at the (a, b) and {b, a) 
position, where it equals one. Since V has at most two nonzero elements, the trace {GVpG 
contains at most 2^N terms. Furthermore, each term is a product of A: + 1 entries of G, e.g. 
Tr GVG = 'P2ki^kbGab + GkaGbk)- We first derive a bound on the resolvent entries Gjk{E + ip) 
down to p > when 9°''^Hg is h-general. Using the delocalization of the eigenvectors we 

have, 


\GjkiE + ip)\ < 


N 


E 


Ui{j)ui{k)\ 
|Ai - z\ 


N 

< GN-^+^ 

i=l 


1 

|Ai - z 
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Define a dyadic decomposition: 


Uo = {j : |A,- -E\< = {j : < \\j - E\} 

Un = {j : < \Xj -E\< 2"iV-^+^}, 1 < n < [loga N]. 


For (5-general 9°'^Hs we have |f7„| < (72”for 0 < n < [log 2 N~\. We can divide the summation 
over i into UnUn, 


\GjkiE + iv)\ <CN-^+^ 




Therefore we have 

F{N-^\Trd[l'>G\ < (7iV3(fc+i)5^ > ^ ^ 

When 6°‘^Hs is not (5-general we still have the deterministic upper bound 

^-ilTra^GI < (7iV3(i+<5)^ 

Therefore we can take E{A) = N~^ Tr(2l — E — ir])~^ for N~^~^ < rj in Lemma 4.3. For B we 
have the upper bound 


B < sup E 

0<s<t,l<i,j<N 
<^(j ]\jl+{3Co+18)5—a 


{N^\hij{s) - f\^ + N\hij{s) - /|)(7iV(3Co + 18)5l ^ (^^3(1+5)^ 


-D 


which yields (5.5). 


□ 


Once we have the above lemma, the following theorem from [15, Theorem 2.1] transforms 
the information of the Green function to the correlation functions of H and will complete the 
proof of Theorem 2.8. 


Theorem 5.3. Let G{z) and Gt{z) denote the Green function of the two matrices H and Eft, 
respectively. Suppose that (5.5) holds for the two random matrices H and Ht for any t = , 

(n) in) 

where e > 0 can be arbitrarily small. Let and be the n-point correlation functions of 
the eigenvalues w.r.t. the probability laws of the matrices H and Ht. Then for any k > 0, 
E S (—2 E k,2 — k), and test function O G G)?°(M”) we have 


lim 

N^OO 



,a, 


.) {E ( 


E + 


CXl 

iV’' 


'E + 


CXr. 


N 


(n) 

Pni 


E + 


Oi\ 

iV 


■■■^'+w)} = 


0 . 


provided that t is chosen so that GtN^^'^^ °‘ < N . 


5.1 Proof of Theorem 2.8 

The universality of the gap statistics follows from Theorem 3.1 and Lemma 5.1. The universality 
of the correlation functions follows from Theorem 3.1, Lemma 5.2 and Theorem 5.3. □ 


6 Universality of sparse random matrices 

In this section we prove Theorem 2.2, the bulk universality of sparse matrices, by checking 
that sparse matrices satisfy the hypotheses of Dehnition 2.4 of stable random matrices. In the 
following we collect some facts about sparse matrices proved in [11]. 

Theorem 6.1. Let H be a sparse Wigner matrix as in Definition 2.1. We denote the eigenvalues 
of H as Xi < X 2 <■■■ < Xn, the corresponding eigenvectors ui,U 2 ,--- ,un, and the resolvent 
G{E -|- ip) = [H — E — ip)~^. Then for all (small) u},5 > 0, (large) D > 0 and large enough 
N > N{oj,5,D) the following holds with probability larger than 1 — N~^. 
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(i) All eigenvalues of H are in the interval [—3, 3], i.e. —3 < Ai < A 2 • • • < Xn-i < 3, except 
for the largest eigenvalue Xn- 

(ii) We have the weak local semi-circle law for individual resolvent entries: there exists a 
constant C such that \Gjk{E + ir])\ < C, for rj > 

(iii) We have the local semi-circle law for the Stieltjes transform of eigenvalues of H: 

\mN{E + ir]) - msciE + ir])\ < N‘^ (- + , \E\ < 5, < p < 10. (6.1) 

\q NrjJ 

Remark 6.2. In fact Theorem 6.1 is still true for generalized sparse matrices, in which we allow 
the variance of each entry K[{hij — /)^] = Sij to be different, given that; (1) they are of the 
same order, i.e. there exists some constants ci and C 2 , such that ciN~^ < Sij < C 2 N~^, (2) 

Theorem 6.3. Sparse random matrices in the sense of Definition 2.1 are stable. Therefore, 
Theorem 2.2 holds. 

Proof. Conditions (1) and (3) of Definition 2.4, i.e., the independent entries and moment con¬ 
ditions, follow from the dehnition of a sparse matrix. The fact that is d-general, i.e., 

condition (4), will follow from the bounds (ii) in Theorem 6.1. If the resolvent elements of 
9°'^Hs are bounded down to the scale rj > then by taking E = Xi and rj = in the 

following identity 


C > Im Gjj{E + ir]) = Y^ (\ 


( 6 . 2 ) 


we see that \ui{j)\‘^ < GN and so the eigenvectors of are completely delocalized. For 

any interval I = \E — ri,E -\- rj\, such that < r/, we have 


N 


C > Im m{E -\-ir]) = — ^ 


> 


1 V 

M 


N^iXi- EY + r/2 - iV (A* - EY + rY - |/|Ar 

2 — 1 Xi^l 




Therefore, we get that ff{i : A, £ /} < G\I\N by rearranging the above expression. 

It therefore suffices to prove that the resolvent entries of the deformed matrix 9°‘^Hs are 
bounded down to the scale rj = For each s, Hg is a sparse matrix and so by Theorem 

6.1 we know that its resolvent entries are bounded with probability greater than 1 — N~^ for 
any large D. 

The deformed matrix 9°‘^Hg is a rank two perturbation of Hg, i.e. 9°‘^Hg = Hg — V, where V 
is the matrix whose matrix elements are zero everywhere except at the (a, b) and (b, a) position, 
where it equals (1 — 9^^)hab(s). By our assumption on the moments of sparse Wigner matrix 
(2.9), we have that 

P (|(1 - 9fYhab(s)j < > 1 - (6.3) 


for any large D. We denote the resolvent of Hg as G = {Hg — z) The resolvent elements of 
9°'^Ha are given by the formula 


W^Hg 



= Gjk + {GV{9^^Hg-z)-Yjk 

<G + 2GN-^G max \{9^^Hg-z)-} 

l<n,l<N 
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with probability greater than 1 — N ^, for N sufficiently large. We get maxj^fc — < 

2C, by taking maximum over j and k in above estimate, and rearranging it. This finishes the 
proof that with probability larger than 1 — N~^ for any (large) number D, is general. 

For the local semi-circle law of i.e., condition (2), the variance of the diagonal terms 

and off-diagonal terms of are e~^N~^ and (1 -|- e~^)j2N respectively. Therefore after 
normalizing by (1 — (1 — e~*))~^ = 1 -|- 0{t), it is a generalized sparse matrix. By (3.12), 
the normalization factor gives us an error of order at most 0{y/i): with probability larger than 
1 - N-^, 


m^^\z) - msc{z) 




\E\ < <r]<W. 


where m^^(z) is the Stieltjes transform of the empirical eigenvalue distribution of There¬ 

fore we can take c((5) = ^ min{a, 6, in (2.19). 

We have checked that sparse random matrices are stable, and so Theorem 2.2 now follows 
from Theorem 2.8. □ 
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